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Abstract. Guided by the relationship between the breadth-first walk of a rooted tree and its 
sequence of generation sizes, we are able to include immigration in the Lamperti representation of 
continuous-state branching processes. We provide a representation of continuous-state branching 
processes with immigration by solving a random ordinary differential equation driven by a pair 
of independent Levy processes. Stability of the solutions is studied and gives, in particular, 
limit theorems (of a type previously studied by Grimvall, Kawazu and Watanabe, and Li) and a 
simulation scheme for continuous-state branching processes with immigration. We further apply 
our stability analysis to extend Pitman's limit theorem concerning Galton- Watson processes 
conditioned on total population size to more general offspring laws. 



1. Introduction 

1.1. Motivation. In this document, we extend the Lamperti representation of Continuous State 
Branching Processes so that it allows immigration. First, we will see how to find discrete (and 
simpler) counterparts to our results in terms of the familiar Galton- Watson process with immi- 
gration and its representation using two independent random walks. 

Consider a genealogical structure with immigration such as the one depicted in Figure 1 . When 
ordering its elements in breadth-first order (with the accounting policy of numbering immigrants 
after the established population in each generation), Xi will denote the number of children of 
individual i. Define a first version of the breadth- first walk x = (xj) by 

xo = and Xi+i = Xi + Xi+i- 

Consider also the immigration process y = (yn)„>o where y„ is the quantity of immigrants arriving 
at generations less than or equal to n (not counting the initial members of the population as 
immigrants). Finally, suppose the initial population has k members. If Cn denotes the number of 
individuals of generations to n, Cn+i is obtained from c„ by adding the quantity of sons of each 
member of the n-th generation plus the immigrants, leading to 

Cn+l = C„ + {Xc„-i+l H K Xc„) + iVn+l " Vn) ■ 

By induction we get 

Cn+l = k + Xc„+ Vn+l- 

Let Zn denote the number of individuals of generation n so that zq = cq = k and for n > 1 

if r/i = Xi — 1 , we can define a second version of the breadth- first walk of the population by setting 

xo = and Xj = Xj_i -|- r]i 
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Figure 1. A genealogical structure allowing immigration 
(so that Xi = Xi — i). We then obtain 

(1) Zn+l = k + Xc„+ Vn+l- 

This representation of the sequence of generation sizes z in terms of the breadth-first walk x and 
the immigration function y can be seen as a discrete Lamperti transformation. It is the discrete 
form of the result we aim at analyzing. However, we wish to consider a random genealogical 
structure which is not discrete. Randomness will be captured by making the quantity of sons 
of individuals an iid sequence independent of the iid sequence of immigrants per generation, 
so that the model corresponds to a Galton- Watson with immigration. Hence x and y would 
become two independent random walks, whose jumps take values in {—1, 0, 1, . . .} and {0, 1, . . .} 
respectively. Discussion of non-discreteness in the random genealogy model would take us far 
apart (we are motivated by Levy trees with or without immigration, discussed for example by 
Duquesne and Le Gall (2002); Lambert (2002); Duquesne (2009); Abraham and Delmas (2009)). 
We only mention that continuum trees are usually defined through a continuum analogue of the 
depth-first walk; our point of view is that generation sizes should be obtained in terms of the 
continuum analogue of the breadth-first walk. Indeed, in analogy with the discrete model, we just 
take X and Y as independent Levy processes, the former without negative jumps (a spectrally 
positive Levy process) and the latter with increasing sample paths (a subordinator). The discrete 
Lamperti transformation of (1) then takes the form 

(2) Z, = x + Xj.^^,^ + Yt. 

This should be the continuum version of a Galton- Watson process with immigration, namely, 
the continuous-state branching processes with immigration introduced by Kawazu and Watanabe 
(1971). 

1.2. Preliminaries. 

1.2.1. (Possibly killed) Levy processes. A spectrally positive Levy process (spLp) is a stochastic 
process X = (Xf)^>Q which starts at zero, takes values on (—00,00], has independent and sta- 
tionary increments, cadlag paths, and no negative jumps. Such a process is characterized by its 
Laplace exponent ^ by means of the formula 

where 

^(X) = -K + aX + + J (^e"^"^ - 1 + Axl^<i j u{dx) 
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for A > 0; here v is the so cahed Levy measure on (0, oo) and satisfies 

1 A v{dx) < oo. 



The constant k will be for us the killing rate; a Levy process with killing rate k can be obtained 
from one with zero killing rate by sending the latter to oo at an independent exponential time of 
parameter n; cr^ is called the diffusion coefficient, while o is the drift. 

We shall also make use of subordinators, which are spLp with increasing trajectories. The 
Laplace exponent $ of a subordinator X is defined as the negative of its Laplace exponent as a 
spLp: 



Since the Levy measure u of a. subordinator actually satisfies 

1 A X u{dx) < oo, 

and subordinators have no Brownian component (i.e. cr^ = 0), we can write 

^'(A) = K + dX + j {l- e^^^') v{dx) . 

So, we have the relationship 

—d = a+ j X v{dx) 
Jo 

between the parameters of X seen as a spLp and as a subordinator. 

1.2.2. Continuous-State Branching processes and the Lamperti representation. Continuous-state 
branching (CB) processes are the continuous time and space version of Galton- Watson processes. 
They were introduced in different levels of generality by Jifina (1958), Lamperti (1967b), and 
Silverstein (1968). They are Feller processes with state-space [0, oo] (with any metric that makes it 
homeomorphic to [0, 1]) satisfying the following branching property: the sum of two independent 
copies started at x and y has the law of the process started at x -|- y. The states and oo are 
absorbing. The branching property can be recast by stating that the logarithm of the Laplace 
transform of the transition semigroup is given by a linear transformation of the initial state. 

As shown by Silverstein (1968), CB processes are in one to one correspondence with Laplace 
exponents of (killed) spectrally positive Levy processes, which are called the branching mecha- 
nisms. In short, the logarithmic derivative of the semigroup of a CB process at zero applied to the 
function x i— )• e~^^ exists and is equal to x i— )• x^'(A). The function ^' is the called the branching 
mechanism of the CB process and it is the Laplace exponent of a spLp. A probabilistic form of 
this assertion is given by Lamperti (1967a) who states that if X is a spLp with Laplace exponent 
^, and for x > we set T for its hitting time of — x, 

= / — — ^7 ds 



JQ X + Xg/s^T 

and C equal to its right-continuous inverse, then 

Zt = X + Xct^j. 

is a CB process with branching mechanism ^, or CB(\I'). This does not seem to be directly 
related to (2). The fact that it is related gives us what we think is the right perspective on the 
Lamperti transformation and the generalization considered in this work. Indeed, as previously 
shown in (Ethier and Kurtz, 1986, Chapter 6§1), Z is the only process satisfying 

(3) Zt = x + Xjt^^^^ 
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which is absorbed at zero . This is (2) in the absence of immigration. To see that a process 
satisfying (3) can be obtained as the Lamperti transform of X, note that if Ct = Jq Zg ds, then 
while Z has not reached zero C is strictly increasing so that it has an inverse, say /, whose 
right-hand derivative is given by 

r ft) = ^— = ± = = ^— 

+ C'4lt) Zj, x + Xcoiit) x + Xt' 

1.2.3. Continuous-State Branching processes with Immigration. Continuous-State Branching Pro- 
cesses with Immigration (or CBI processes) are the continuous time and space version of Galton- 
Watson processes with immigration and were introduced by Kawazu and Watanabe (1971). They 
are Feller processes with state-space [0, oo] such that the logarithm of the Laplace of the transition 
semigroup is given by an affine transformation of the initial state. (They thus form part of the 
affine processes studied by Dawson and Li (2006).) As shown by Kawazu and Watanabe (1971), 
they are characterized by the Laplace exponents of a spLp and of a subordinator: the logarithmic 
derivative of the semigroup of a CB process at zero applied to the function x i— )• e~^^ exists and 
is equal to the function 

X x^{X) - $(A) 

where ^ is the Laplace exponent of a spLp and ^ is the Laplace exponent of a subordinator. 
They are respectively called the branching and immigration mechanisms and characterize the 
process which is therefore named CBI(^, $). 

We aim at a probabilistic representation of CBI processes in the spirit of the Lamperti repre- 
sentation. 

1.3. Statement of the results. We propose to construct a CB^^*, <&) that starts at x by solving 
the functional equation 

(4) Z, = x + Xj.^^,^ + Yt. 

We call such a process Z the Lamperti transform of {X, x + Y) and denote it by Z = L{X, x + Y); 
however, the first thing to do is to show that there exists a unique process which satisfies (4). 
When Y is zero, a particular solution to (4) is the Lamperti transform of X + x recalled above. 
Even in this case there could be many solutions to (4), in clear contrast to the discrete case 
where one can proceed recursively to construct the unique solution. Our stepping stone for the 
general analysis of (4) is the following partial result concerning existence and uniqueness proved 
in Section 2. 

A pair of cadlag functions (/, g) such that / has no negative jumps, g is non-decreasing and 
/(O) +g(0) > is termed an admissible breadth-first pair; f and g will be termed the reproduction 
and immigration functions respectively. When g is constant, we say that f -\-g is absorbed at zero 
if f{x) + (7 = implies f{y) + 5 = for all y > x. 

Theorem 1. Let {f,g) be an admissible breadth- first pair. There exists a nonnegative h satisfying 
the equation 

Kt) = f( [ h{s) ds)+git). 



\Jo / 

Furthermore, the solution is unique when g is strictly increasing, when f + (^(0) is a strictly 
positive function, or when g is constant and f + g is absorbed at zero. 

In the context of Theorem 1, much is gained by introducing the function c given by 

c{t) = I h{s) ds, 
Jo 
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which has a right-hand derivative equal to h. This is because the functional equation for h 
can then be recast as the initial value problem 



IVP(/,<7) 



c'+ = f OC + g 
c(0) = 



Our forthcoming approximation results for the function h of Theorem 1 rely on the study of 
a functional inequality. Let (/, g) be an admissible breadth-first pair. We will be interested in 
functions c which satisfy 

(5) / /_ o c(r) + g{r) dr < c{t) — c{s) < f o c(r) + g{r) dr, for s < t. 

J s J s 

Note that any solution c to IVP(/, g) satisfies (5): the second inequality is actually an equality by 
definition of IVP(/, (7) and since f > f- as f has no negative jumps, we get the first inequality. 
Hence, the functional inequality (5) admits solutions. Regarding uniqueness, if the solution to 

(5) is unique, then the solution to IVP(/, is unique and since the latter is non-negative and 
non-decreasing, so is the former. Also, similar sufficient conditions for uniqueness of IVP(/, (7) of 
Theorem 1 imply uniqueness of non-decreasing solutions of the functional inequality (5). 

Proposition 1. Let {f,g) be an admissible breadth-first pair. If either g is strictly increasing, 
/_ + g{0) is strictly positive, or g is constant and /_ + g{0) is absorbed at zero, then (5) has an 
unique non- decreasing solution starting at zero. 

However, as is shown in Subsection 4.1, assuming that (5) admits an unique solution is stronger 
than just assuming that IVP(/, 5) has an unique solution. 

As a consequence of the analytic Theorem 1, we solve a probabilistic question raised by Lambert 
(1999, 2007). 

Corollary 1. Let X be a spectrally positive a-stable Levy process. For any cadlag and strictly 
increasing process Y independent of X , there is weak existence and uniqueness for the stochastic 
differential equation 

(6) Zt = x+ f \z,\ V"dx, + yt. 



When X is twice a Brownian motion and Yt = 5t for some 5 > this might be one of the 
simplest proofs available of weak existence and uniqueness of the SDE defining squared Bessel 
processes, since it makes no mention of the Tanaka formula or local times; it is based on Knight's 
theorem and Theorem 1. When X is a Brownian motion and dYt = b{t) dt for some Lipschitz and 
deterministic b : [0, 00) — )• [0, 00), Le Gall (1983) actually proves pathwise uniqueness through a 
local time argument. Our result further shows that if b is measurable and strictly positive then 
there is weak uniqueness. In the case Y is an (a — l)-stable subordinator independent of X, we 
quote Lambert (1999, 2007) 

..whether or not uniqueness holds for (6) remains an open question. 

Corollary 1 answers affirmatively. Note that when Y = the stated result follows from Zanzotto 
(2002), and is handled by a time-change akin to the Lamperti transformation. Fu and Li (2010) 
obtain strong existence and pathwise uniqueness for a different kind of SDE related to CBI 
processes with stable reproduction and immigration. 

Regarding solutions to (4), Theorem 1 is enough to obtain the process Z when the subordinator 

Y is strictly increasing. When Y is compound Poisson, a solution to (4) can be obtained by pasting 
together Lamperti transforms. However, further analysis using the pathwise behavior of X when 

Y is zero or compound Poisson implies the following result. 
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Proposition 2. Let x > 0, X be a spectrally positive Levy process and Y an independent subor- 
dinator. Then there is a unique cddldg process Z which satisfies 

The above equation is satisfied by any cddldg process Z satisfying the functional inequality 

x + Xj^^^^^_+Yt<Zt<x + Xj. + Yt 

which also has a unique solution. 

Our main result, a pathwise construction of a CBI(^, $), is the following. 

Theorem 2. Let X be a spectrally positive Levy process with Laplace exponent ^ and Y an 
independent subordinator with Laplace exponent The unique stochastic process Z which solves 



Zt = x + Xjt^^^^ + Yt 



is a CBI(^, ^) that starts at x. 



We view Theorems 1 and 2 as a first step in the construction of branching processes with 
immigration where the immigration can depend on the current value of the population. One 
generalization would be to consider solutions to 

Zt=X + Xjt ^^^2^-, + Yjt y^^^z,)ds^ 

where a is interpreted as the breeding rate and h as the rate at which the arriving immigration is 
incorporated into the population. For example Abraham and Delmas (2009) consider a continuous 
branching process where immigration is proportional to the current state of the population. This 
could be modeled by the equation 

which, thanks to the particular case of Theorem 2 stated by Lamperti (1967a), has the law 
of a CB(^ — a$) started at x; this is the conclusion of Abraham and Delmas (2009), where 
they rigorously define the model in terms of a Poissonian construction of a more general class 
of CBI processes which is inspired in previous work of Pitman and Yor (1982) for CBIs with 
continuous sample paths. Another representation of CBI processes, this time in terms of solutions 
to stochastic differential equations was given by Dawson and Li (2006) under moment conditions. 

The usefulness of Theorem 2 is two-fold: firstly, we can use known sample path properties of X 
and Y to deduce sample-path properties of Z, and secondly, this representation gives a particular 
coupling with monotonicity properties which are useful in limit theorems involving Z, as seen in 
Corollary 6, Corollary 7, and Theorem 4. Simple applications of Theorem 2 include the following. 

Corollary 2 (Kawazu and Watanabe (1971)). If^ is the Laplace exponent of a spectrally positive 
Levy process and $ is the Laplace exponent of a subordinator, there exists a CBI process with 
branching mechanism ^ and immigration mechanism <I>. 

Corollary 3. A CBI(^',<^) process does not jump downwards. 

Caballero et al. (2009) give a direct proof of this when <I> = 0. 

Corollary 4. Let Z be a CBI('I',<I>) that starts at x > 0, let <I> be the right- continuous inverse of 
^ and define 

/,x log I log 1 1 

a{t) = . 

^{t-^ log |logt|) 
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There exists a constant C, (in general non-zero) such that 

liminf —. — — = C. 

t^o a{xt) 

The case x = in Corollary 4 is probably very different, as seen when ^'(A) = 2A^ and 
$(A) = dX which corresponds to the squared Bessel process of dimension d. Indeed, (Ito and 
McKean, 1974, p. 80) show that for a squared Bessel process Z of integer dimension that starts 
at 0, we have 

lim sup — r = 1 . 

t^o 2ilog|logt| 

We have not been able to obtain this result using the Lamperti transformation. However, note 
that starting from positive states we can obtain the lower growth rate, since it is the reproduction 
function X that determines it, while starting from 0, it is probably a combination of the local 
growth of X and Y that drives that of Z. 

A solution c to IVP(/, (7) is said to explode if there exists t G (0,oo) such that c(t) = 00. 
(Demographic) explosion is an unavoidable phenomena of IVP(/, (7). When / > and g = 0, it 
is known that explosion occurs if and only if 

dx < 00. 



Actually, even when there is immigration, the main responsible for explosion is the reproduction 
function. 

Proposition 3. Let {f,g) be an admissible pair and let = max(/, 0). 

(1) If f°° l//+(x) dx = 00 then no solution to IVP(/,g) explodes. 

(2) If J l/f^{x) dx < 00, lim^_s.oo /(a;) = 00 and g{oo) exceeds the maximum of —f then 
any solution to TVP(f,g) explodes. 



We call / an explosive reproduction function if 

1 



fix < 00. 



f+{x] 

Recall that 00 is an absorbing state for CBI processes; Proposition 3 has immediate implications 
on how a CBI process might reach it. First of all, CBI processes might jump to 00, which happens 
if and only if either the branching or the immigration correspond to killed Levy processes. When 
there is no immigration and the branching mechanism ^ has no killing rate, the criterion is due 
to Ogura (1970) and Grey (1974), who assert that the probability that a CB(\I') started from 
X > is absorbed at infinity in finite time is positive if and only if 

dX > —00. 



M/(A) 

One can even obtain a formula for the distribution of its explosion time (cf. the proof of Theorem 
2.2.3.2 in (Lambert, 2008, p. 95)). We call such an explosive branching mechanism. From 
Proposition 3 and Theorem 2 we get: 

Corollary 5. Let x > 0. 

(1) The probability that a CBI(^, $) Z that starts at x jumps to 00 is positive if and only if 
^'(0) or <I*(0) are non-zero. 

(2) The probability that Z reaches 00 continuously is positive if and only if ^(0) = and "if 
is an explosive branching mechanism. 
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(3) The probability that Z reaches oo continuously is equal to 1 if^{0) = $(0) =0, ^ is not 
zero and ^ is explosive. 

We mainly use stochastic integration by parts in our proof of Theorem 2; however, a weak 
convergence type of proof, fohowing the case <I> = presented by Cabahero et al. (2009), could 
also be achieved in conjunction with a stability result, based on the forthcoming Theorem 3. 

The following result deals with stability of IVP(/, g) under changes in / and g and even includ- 
ing a discretization of the initial value problem itself. Indeed, consider the following approximation 
procedure: given cr > 0, called the span, consider the partition 

ti = ia, i = 0,1,2, ... , 

and construct a function c'^ by the recursion 

c'^(O) = 

and for t £ [ti-i,ti): 

c^t) = c^ti.i) + {t- [/ o c^(t,„i) + <7(ti_i)]+ . 
Equivalently, the function c"^ is the unique solution to the equation 

IVP,(/,5) iC^W = f[foc'^{[s/a\a) + g{[s/a\a)]-^ ds. 
Jo 

We will write IVPo(/, 5) to mean TVP{f,g). Let 1)+ denote the right-hand derivative. 

The stability result is stated in terms of the usual Skorohod Ji topology for cadlag functions: a 
sequence converges to / if there exist a sequence of homeomorphisms of [0, 00) into itself such 
that 

fn — f ° \i and A„ — Id converge to zero uniformly on compact sets 

(where Id denotes the identity function on [0,oo)). However, part of the theorem uses another 
topology on non-negative cadlag functions introduced by Caballero et al. (2009), which we propose 
to call the uniform Ji topology. Consider a distance d on [0, 00] which makes it homeomorphic 
to [0, 1]. Then the uniform Ji topology is characterized by: a sequence fn converges to / if there 
exist a sequence of homeomorphisms of [0, 00) into itself such that 

d-ifn, f ° An) — 5- and A„ — Id — )• uniformly on [0, 00). 

Theorem 3. Let (/, g) be an admissible breadth-first pair and suppose there is a unique non- 
decreasing function c which satisfies c(0) = and (5) (and is therefore the unique solution to 
IVP{f,g)); define its explosion time by 

r = inf {t > : c{t) = 00} G (0,oo]. 

Let {fn,gn) be admissible breadth-first pairs. Suppose fn^f o,nd gn ^ g in the Skorohod Ji 
topology and that Un is a sequence of non-negarive real numbers which tend to zero. Let Cn be 
the unique solution to IVPa„ifn, dn) when cr„ > and any solution to WP{fn, dn) when an = 0. 
Then Cn ^ c pointwise and uniformly on compact sets of [0,r). 

Furthermore, if f o c and g do not jump at the same time then -D+c^ — t- D+c 

(1) in the Skorohod Ji topology if t = 00, and 

(2) in the uniform Ji topology i/r < 00 if we additionally assume that fn{x) , f{x) — )• 00 as 

— )• 00 uniformly in n. 

It is not very hard to show that the jumping condition of Theorem 3 holds in a stochastic 
setting. 
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Proposition 4. Let X he a spLp , Y an independent subordinator with Laplace exponents ^ and 
$ and, for x >0, let Z the unique process such that 

Zt = X + Xct + Yt where Ct = Zg ds. 

Jo 

Almost surely, the processes X o C and Y do not jump at the same time. 

From Theorem 3 and Propositions 2 and 4, we deduce tlie following weak continuity result. 

Corollary 6. Let ^n-,^ be Laplace exponents of spLps and ^n,^ be Laplace exponents of sub- 
ordinators and suppose that — ^ ^ o-n-d ^ ^ pointwise. If is a sequence in [0, oo] 
converging to x and Zn (resp. Z) are CBIs with branching and immigration mechanisms o-nd 
$n (resp. ^ and ^) and starting at x„, (resp. x) then Zn ^ Z in the Skorohod Ji topology on 
cadldg paths on [0, oo] if is non-explosive and in the uniform Ji topology if ^ is explosive. 

Theorem 3 also allows us to simulate CBI processes. Indeed, if we can simulate random 
variables with distribution Xt and Yt for every t > 0, we can then approximately simulate the 
process Z as the right-hand derivative of the solution to IVPo-(X, x + Y). (Alternatively, if we can 
approximate X and Y, for example by compound Poisson processes with drift, we can also apply 
IVPo- to approximate the paths of Z.) The procedure IVPo-(^, x + Y) actually corresponds to an 
Euler method of span a to solve IVP(X, x + Y). Theorem 3 implies the convergence of the Euler 
method as the span goes to zero when applied to IVF{X,x + Y), even with the discontinuous 
driving functions X and Y\ 

We also give an application of Theorem 3 to limits of Galton- Watson processes with immigra- 
tion. Let X" and Y"" are independent random walks with step distributions //„ and Vn supported 
on{— 1,0,1,...} and {0, 1, 2, . . .} and for any /c„ > 0, define recursively the sequences and 
by setting 

As discussed in Subsection 1.1, the sequence is a Galton- Watson process with immigration 
with offspring and immigration distributions /x„ and f„. However, if X" and Y^ are extended by 
constancy on [m,m + 1) for m > (keeping the same notation), then is the approximation 
of the Lamperti transformation with span 1 applied to X"' and Y"" and Z" is the right-hand 
derivative of C". In order to apply Theorem 3 to these processes, define the scaling operators 
by 

S'afit) = If {at) . 



Corollary 7. Suppose the existence of sequences an,bn such that 

X" /n and Y^" In 

converge weakly to the infinitely divisible distributions /x and v corresponding to a spectrally pos- 
itive Levy process and a subordinator; denote by ^ and $ their Laplace exponents. Suppose that 
6„ — )■ oo and, for any a > 0, a^an\/n' ~^ oo. Let /c„ — )• oo and suppose that either 

— — — )• c G [0, oo) or - — — — > 

as n ^ oo. Setting = in the first case and Cn = xa^f^^/^^/kn in the second, we have 

that 
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converges in distribution, towards a CBI(c^', in the first case and towards a CB(^') in the 
second. The convergence takes place in the Skorohod Ji topology if^ is non-explosive and in the 
uniform Ji topology otherwise. 

When ^ is non-explosive and $ = 0, the above theorem was proved by Grimvall (1974). He 
also proved the convergence of finite-dimensional distributions in the explosive case, which we 
complement with a limit theorem. For general but non-explosive ^, a similar result was proven 
by Li (2006). However, as will be seen in the proof (which relies on the stability of the Lamperti 
transformation stated in Theorem 3), if the convergence of S^^X^ and S^^Y^ takes place almost 
surely, then S^^^Z'^ also converges almost surely. 

The stability result of Theorem 3 applies not only in the Markovian case of CBI processes. As 
an example, we generalize work of Pitman (1999) who considers the scaling limits of conditioned 
Galton- Watson processes in the case of the Poisson offspring distribution. Let ^ be an offspring 
distribution with mean 1 and suppose that Z"' is a Galton- Watson process started at kn and 
conditioned on 



We shall consider the scaling limit of as n — )• cx) whenever the shifted reproduction law 
Afc = IJ-k+i is in the domain of attraction of a stable law without the need of centering. The 
scaling limit of a random walk with step distribution jl is then a spectrally positive stable law of 
index a G (1, 2] with which one can define, for every I > the first passage bridge starting at 
I and ending at of length 1 of the associated Levy process. Informally this is the stable process 
started at I, conditioned to be above on [0, 1] and conditioned to end at at time 1. This 
intuitive notion was formalized by Chaumont and Pardo (2009). The Lamperti transform of 
will be the right hand derivative of the unique solution to IVP(F',0). 

Theorem 4. Let be a Galton-Watson process with critical offspring law /i which starts at kn 
and is conditioned on ^"^i Zf = n. Let S be a random walk with step distribution ^ and suppose 
there exist constants a„ — )• oo such that (Sn — n) /a„ converges in law to a spectrally positive 
stable distribution with Laplace exponent ^ . Let X be a Levy process with Laplace exponent ^ 
and F^ its first passage bridge from I > to of length 1. If kn/a^ — )• / then the sequence 

n/a„ 

converges in law to the Lamperti transform of F^ in the Skorohod Ji topology. 

When a = 2, the process F^ is a Bessel bridge of dimension 3 between / and of length 1, 
up to a normalization factor. In this case, (Pitman, 1999, Lemma 14) tells us that the Lamperti 
transform Z^ of F^ satisfies the SDE 



driven by a Brownian motion and it is through stability theory for SDEs that Pitman (1999) 
obtains Theorem 4 when /i is a Poisson distribution with mean 1. Theorem 4 is a complement 
to the convergence of Galton-Watson forests conditioned on their total size and number of trees 
given in Chaumont and Pardo (2009). When / = 0, our techniques cease to work. Indeed, 
the corresponding process F^ would be a normalized Brownian excursion above zero and the 
problem IVP(F'',0) does not have a unique solution, as discussed at the beginning of Section 
2. Hence, even if our techniques yield tightness in the corresponding limit theorem with I = 0, 
we would have to give further arguments to prove that any subsequential limit is the correct 



oo 
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solution IVP(-F'^, O). The limit theorem when / = and a = 2 was conjectured by Aldous (1991), 
and proved by Drmota and Gittenberger (1997) by analytic methods. For any a G (1,2], the 
corresponding statement was stated and proved by Kersting (1998) by working with the usual 
Lamperti transformation, which chooses a particular solution to IVP(i^'',0). 

The paper is organized as follows. Theorem 1, Proposition 2, and Corollary 1 are proved in 
Section 2 which focuses on the analytic aspects of the Lamperti transformation and its basic 
probabilistic implications. The representation CBI processes of Theorem 2 is then proved in 
Section 3, together with Proposition 4, Corollary 4 and Corollary 5. Finally, Section 4 is devoted 
to the stability of the Lamperti transformation with a proof of Theorem 3, Proposition 1, Corollary 
6, Corollary 7, and Theorem 4. (Corollaries 2 and 3 are considered to follow immediately from 
Theorem 2; proofs have been omitted.) 

2. The generalized Lamperti transformation as an initial value problem 

Let (/, g) be an admissible breadth-first pair, meaning that / and g are cadlag functions with 
g increasing, / without negative jumps and /(O) + gi^) > 0. We begin by studying the existence 
of a nonnegative cadlag function h which satisfies 

(7) h{t) = f([\{s)ds]+g{t); 



\Jo / 

a priori there might be many solutions. 

When g is identically equal to zero, a solution is found by the method of time-changes: let r 
be the first hitting time of zero by /, let 

^* = / ttA — V 
Jo /(sAr) 

and consider its right-continuous inverse c so that 

h = f oc 

satisfies (7) with g = and it is the only solution for which zero is absorbing. A generalization of 
this argument is found in (Ethier and Kurtz, 1986, Chapter 6§1). In this case the transformation 
which takes f to h is called the Lamperti transformation, introduced by Lamperti (1967a). There 
is a slight catch: if / is never zero and goes to infinity, then h exists up to a given time (which 
might be infinite) when it also goes to infinity. After this time, which we call the explosion time, 
we set h = oo. With this definition, note that c and h become infinite at the same time. 

Solutions to (7) are not unique even when g = as the next example shows: take f{x) = 
Y^|l — x| , / > 0, and consider 



(2 - t)+ 
hi{t) = - — and h2{t) 



^ if t < 2 
if 2 < t < 2 / 
iit>2 + l 



2 

Then hi and /12 are both solutions to (7). As discussed in the introduction, a probabilistically 
relevant example of non-uniqueness is obtained when g = and / is the typical sample path of 
a normalized brownian excursion e = {et,t > 1). (See Chapter 11§3 of Revuz and Yor (1999) for 
its definition as a 3-dimensional Bessel bridge.) Indeed, with probability 1, e has a continuous 
trajectory which is positive exactly on (0, 1). Hence, is a solution to IVP(e, 0). However, its link 
with the 3-dimensional Bessel process (and time reversal) allows one to prove that ^/s = o{es) as 
s — >■ 0-|- (and a corresponding statement as s — 1— ) so that almost surely 

— as < 00. 
e,9 
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Hence, one can define the Lamperti transform of e, which is a non-trivial solution to IVP(e,0). 
The Lamperti transformation is well defined under more general excursion laws as discussed by 
Miermont (2003). 

We propose to prove Theorem 1 by the following method: we first use the solution for the 
case 5 = to establish the theorem when g is piecewise constant. When g is strictly increasing, 
we approximate it by a strictly decreasing sequence of piecewise constant functions g-n > 9 and 
let hn be the solution to (7) which uses We then consider the primitive c„ of /i„ starting at 
zero, show that it converges, and this is enough to prove the existence of a function whose right- 
continuous derivative exists and solves (7). Actually, it is using primitives that one can compare 
the different solutions to (7) (and study uniqueness) and this is the point of view adopted in what 
follows. To this end, we generalize (7) into an initial value problem for the function c: 



IVP(/,5,x) 



c'+{t) = f o dt) + g{t) 
c(0) = X 



(The most important case for us is x = and we will write IVP(/, 5) when referring to it.) We 
shall term 

• / the reproduction function, 

• g the immigration function, 

• X the initial cumulative population, and 

• c will be called the cumulative population. 

• c'^ is called the population profile. 

• A solution c to TVP{f,g,x) is said to have no spontaneous generation if the condition 
c'_|.(i) = implies that c{t + s) = c{t) as long as g{t + s) = g{t). 

In the setting of Theorem 1, spontaneous generation is only relevant when g is piecewise constant 
and it will be the guiding principle to chose solutions in this case. 

A solution to TVP{f,g,x) without spontaneous generation when g' is a constant 7 is obtained 
by setting fx{s) = f{x + s) + 7, calling the Lamperti transform of fx and setting 

ft 



We then have: 



ct = X + / hx{s) ds. 
Jo 

c+{t) = hx{t) = fx(^j\x{s) ds^ =f[x + j\x{s) ds^ + j = f{c{t))+g{t). 
Let g be piecewise constant, say 

n 

1=1 

with 71 < 72 < • • • < 7n and = to < ^i < ■ ■ " < ^n- Let us solve (7) by pasting the solutions on 
each interval: let ^1 solve IVP(/, 71,0) on [0,ti] without spontaneous generation. Let c equal ipi 
on [0,ti]. Now, let V2 solve IVP(/, 72, c(ii)) without spontaneous generation. (If c(ti) = 00, we 
set 11)2 = 00.) Set c{t) = ip2{t — ti) for t G [^1,^2] so that c is continuous. Also, for t G [^1,^2] we 
have 

cV(t) = M - ti) = fiMt - h)) + 72 = /(c(t)) + 9it) . 
We continue in this manner. Note that if c'_^_ reaches zero in [ti-i,ti), say at t then c is constant 
on [t,ti) and that solves (7) when g is piecewise constant. By uniqueness of solutions to (7) 
which are absorbing at zero when g = 0, we deduce the uniqueness of solutions to IVP(/, 5, 0) 
without spontaneous generation when the immigration is piecewise constant. 
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We first tackle the non-negativity assertion of Theorem 1. Since / is only defined on [0,oo), 
negative values of c do not make sense in equation (7). One possible solution is to extend / to M 
by setting f{x) = /(O) for x < 0. 

Lemma 1. Any solution h to (7) is non-negative. 

Proof. Let h solve (7) where / is extended by constancy on (— oo,0] and define 

c{t) = / h{s) ds, 
Jo 

so that c solves IVP{f,g). We prove that /i > by contradiction. Assume there exists t > 
such that h{t) < 0. Note that since h has no negative jumps, h can only reach negative values 
continuously and, since h is right-continuous, if it is negative at a given t, then there exists t' > t 
such that h is negative on [t,t'). Hence there exists e > such that 

{t>0: h{t) = and /i< on {t, t + e)} ^ 0. 

Let r be its infimum. We assert that r > and c(r) > 0. Indeed, if r = then c would be 
strictly decreasing and negative on (0, e), which would imply that 

h{t) = f o c{t) + g{t) = /(O) + g{t) > /(O) + g{0) > for t G (0, e), 

a contradiction. A similar argument tells us that c(r) > 0. We finish the proof by showing 
the existence of ti < r and t2 G {t,t + e) such that h{ti) > and c(ti) = 0(^2), implying the 
contradiction 

< h{ti) = f o c{ti) + g{ti) = f o c{t2) + g{ti) < f o c{t2) + g{t2) < 0. 

Indeed, given that c(r) > we can assume that c(r + e) > by choosing a smaller e, and then 
let Tl be the last time before r that c is below c(t -|- e) and T2 the first instant after ri that c 
equals c(r). Note that T2 < r. Since 

T2 

h{r) dr = c{t2) — c(ri) = c(r) — c(r -|- e) > 0, 

1 

there exists r G (ri, T2) such that h{r) > and by construction c(r) G c((r, r -|- e)). □ 

2.1. Monotonicity and existence. We now establish a basic comparison lemma for solutions 
to IVP(/, (7) which will lead to the existence assertion of Theorem 1. 

Lemma 2. Let c and c solve IVP(/, (7) and TVF{f,g). If 

5(0) + /(O) < 5(0) + /(O) , /</, g<~g 

and either g^ < g^ or /_ < /_ then ct < Ct for every t that is strictly positive and strictly smaller 
than the explosion time of c. 

It is important to note that the inequality c < c cannot be obtained from the hypothesis 
g < g using the same reproduction function /. Indeed, we would otherwise have uniqueness for 
IVP(/, 5) which, as we have seen, is not the case even when g = Also, since both c and c 
begin at and equal 00 after their explosion time, we always have the inequality c < c under the 
conditions of Lemma 2. 

Proof. Let r = inf {t > : c{t) = c{t)]. Since 

cV(o) = /(o) + 5(0) < /(o) + m = c;(o) , 

and the right-hand derivatives of c and c are right-continuous, then r > and c < c on (0,r). 
Note then that the explosion time of c cannot be smaller than r, since this would force c to 
explode before r and so c would equal c before r. 
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We now argue by contradiction. If r were finite, we know that 

c(r) = c(r) 

leaving us with two cases: 

c(r) = c(r) = oo and c(r) = c(r) < oo. 

In the former, we see that r is the explosion time of c and so the statement of Lemma 2 holds. 
In the latter case, 

c'_(r) = /(c(t) -) + gir-) = f{c{T) -) + g{T-) < /(c(r) -) + g{T-) = (r) . 

It follows that c'_ < c'_ in some interval (r — e, r). However, for < t < r we have c(t) < c(t) 
and this implies the contradiction 

c(r) < c(r) . 

□ 

Proof of Theorem 1, Existence. Consider a sequence of piecewise constant cadlag functions 
satisfying 5n+i(0) < 5n(0), S-n+i- < 5n- and such that gn ^ g pointwise. Let c„ solve IVP(/,5n) 
with no spontaneous generation. By Lemma 2, the sequence of non-negative functions c„ is 
decreasing, so that it converges to a limit c. Let 

r = inf {t > : c{t) = oo} = lim inf {t > : c„(t) = oo} . 

n— >oo 

Since / is right-continuous and c < Cn, f o Cn + gn converges pointwise to f o c + g on [0,t). By 
bounded convergence, for t G [0,t): 

c{t) = lim Cn{t) = lim / / o c„(s) + ^^(s) = / f o c{s) + g{s) ds. 

Hence, h = c'^ proves the existence part of Theorem 1. □ 

2.2. Uniqueness. To study uniqueness of IVP(/, g), we use the following lemma. 

Lemma 3. If g is strictly increasing and c solves TVP(f,g) then c is strictly increasing. 

Proof. Note that by Lemma 1, the right-hand derivative of c is non- negative, so that c is non- 
negative and non-decreasing. By contradiction: if c had an interval of constancy [s, t], with t > s, 
then 

> foc{s) + g{s) 
= 0. 

□ 

Remark. As we shall see in the proof of the uniqueness assertion of Theorem 1 , if we can guar- 
antee that all solutions to IVP{f,g) are strictly increasing then uniqueness holds for IVP(/, 5). 
Note that if / + 5(0) is strictly positive then /(x) + g{t) > for all x > and i > 0, so that all 
solutions to IVP(/, are strictly increasing. 
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Proof of Theorem 1, Uniqueness. Let c and c solve YVF{f,g). To show that c = c, we argue by 
contradiction by studying their inverses i and i. 

Suppose that c and c are strictly increasing. Then i and i are continuous. If c 7^ c then i ^ i 
and we might without loss of generality suppose there is xi such that i{xi) < i{xi). Let 

xq = sup |x < xi : i(x) > 

and note that, by continuity of i and i, xq < x\ and i < i on (xq, xi]. Since i and i are continuous, 
they satisfy 

rv 1 

Jo /(a^) + 9° 

There must exist x £ [xo,xi] such that i'{x) and i'{x) both exist and the former is strictly smaller 
since otherwise the inequality i < i would hold on [xo,a;i]. For this value of x: 

1-1 
fix) = -7- -g° < -TfT-: -go i{x) = f{x) , 
i'{x) v[x) 

which is a contradiction. 

Note that all solutions to IVP(/, (7) are strictly increasing whenever g is strictly increasing (by 
Lemma 3) or / is strictly positive, which implies uniqueness to IVP(/, (7) in these cases. 

When g is constant and f+g is absorbed at 0, meaning that if f{s)+g{Q) = then f{t)+g{0) = 
for all t > s, we can directly use the Lamperti transformation to obtain uniqueness. Indeed, 
solutions to TVP{f,g) do not have spontaneous generation and, as stated in the introduction to 
Section 2 (cf. page 11), there is an unique solution to IVP(/ + 5(0) , 0) without spontaneous 
generation. □ 

2.3. Uniqueness in the stochastic setting. We now verify that solutions to (4) are unique 
even if the subordinator Y is compound Poisson. 

Proof of Proposition 2. Let X be a spLp and Y and independent subordinator. We first prove 
that there is an unique process Z which satisfies 

Zt = x + x(^J^ Z,ds^ + Yt. 

When Y is an infinite activity subordinator (its Levy measure is infinite or equivalently it has 
jumps in any nonempty open interval) or it has positive drift, then its trajectories are strictly 
increasing and so uniqueness holds thanks to Theorem 1. 

It then suffices to consider the case when y is a compound Poisson process. There is a simple 
case we can establish: if X is also a subordinator and x > then all solutions to TVP{X, x + Y) 
are strictly increasing and so uniqueness holds (again by Theorem 1). It remains to consider two 
cases: when X is a subordinator and x = and when X is not a subordinator. In the first, note 
that zero solves IVP(X, 0) and since every solution is non-negative, zero is the smallest one. To 
prove uniqueness, let be the (unique) solution to IVP(X, x), so that is greater than any 
solution to IVP(X, 0) by Lemma 2. If we prove that as x — t- 0, — >• 0, then all solutions to 
IVP{X, 0) are zero and so uniqueness holds. For this, use the fact that as t — >• 0, X^/t converges 
almost to the drift coefficient of X, say d G [0,oo) (cf. (Bertoin, 1996, Ch. Ill, Proposition 8, 
p. 84)) so that 

— ds = 00. 

0+ 
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Let be the (continuous) inverse of (note that is strictly increasing). Since 

we see, by Fatou's lemma, that — )• oo as x — )■ 0, so that — )• 0. Now with X still a 
subordinator and Y compound Poisson, the preceding case implies that the solution to IVF{X, Y) 
is unique until the first jump time of Y; after this jump time, all solutions are strictly increasing, 
hence uniqueness holds. 

The only remaining case is when Y is compound Poisson and X is not a subordinator. The last 
hypothesis implies that is regular for (— oo,0), meaning that on every interval [0, e), X visits 
(-00,0) (cf. (Bertoin, 1996, Ch. VII, Thm. 1, p. 189)). From this, it follows that if T is any 
stopping time with respect to the filtration a{Xs,s < t) V (7{Y) ,t >0, then X visits (— oo,^^) 
on any interval to the right of T. Let C be any solution to IVP(X, x + Y); we will show that it 
has no spontaneous generation. Since there is an unique solution without spontaneous generation 
when Y is piecewise constant (as discussed in the introduction to Section 2), we get uniqueness. 
Let 

[Ti.i,Ti),i = 1,2,... 

be the intervals of constancy of Y; if C has spontaneous generation on one of these, say [Tj_i, Tj), 
then X reaches the level — lTi_i and then increases, which we know does not happen since the 
hitting time of {— lTi_i} by the process X is a stopping time with respect to the filtration 
a{Xs,s < t)\Ja{Y),t> 0. 

We end the proof by showing that any cadlag process Z satisfying 



(8) x + X_[ I Z,ds] +Yt<Zt<x + X( I Zsds] +Yt 

actually satisfies 



Let 



ft 

x + X{ I Zsds] + Yt. 

r-t 

Ct = I Zgds. 



When Y is strictly increasing, an argument similar to the proof of the Monotonicity Lemma 
(Lemma 2) tells us that C is strictly increasing, so that C actually satisfies IVP(X, x + Y). 

When Y = 0, the previous argument shows that, as long as Z has not reached 0, C coincides 
with the solution to TVF{X,x). If Z is such that 

inf {t>0:Zt = 0} = mf{t>0: Zt- = 0} , 

then C solves IVP(X, x), which has an unique solution, so that (8) has an unique solution. We 
then see that the only way in which Z can cease to solve IVP(X, x) is if X is such that 

ro+ = inf {t > : x + = 0} < inf {t > : X + = 0} = Tq, 

which is ruled out almost surely by quasi left-continuity of X. Indeed, ro+ is the increasing limit 
of the stopping times 

Te = '\ni{t>Q:x + Xt<e] 

which satisfy < T^' if e < e' since X has no negative jumps. Hence X is almost surely 
continuous at Tqj^ which says that x + Xtqj^ = almost surely. In the remaining case when Y 
is a (non-zero) compound Poisson process, we condition on Y and argue similarly on constancy 
intervals of y. □ 
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2.4. Explosion. We now turn to the explosion criteria of solutions of IVP(/, of Proposition 
3. 

Proof of Proposition 3 . 

(1) If = oo, let c be any solution to IVP(/, g). We show that c is finite at 

every t > 0. Indeed, using the arguments of Lemma 2, we see that c is bounded by any 
solution to IVP(/, 1 + g{t)) on the interval [0, t]. A particular solution to IVP(/, 1 + g(t)) 
is obtained by taking the right-continuous inverse of 



Since 



dx = oo, 



/o /+(x) + l + <7(t) 
the particular solution we have considered is everywhere finite. 
(2) Let c be a solution to IVP(/, g) where / is an explosive reproduction function, lim^j^oo fix) 
oo and g{oo) exceeds the maximum of — /. To prove that c explodes, choose T > such 
that fix) + g{t) > for all x > and t > T. Then f o c + g > on [T,oo). Let 
M = c{T). We then consider the right-continuous inverse i of c (which is actually an 
inverse on [M, oo)) and note that for y > M 

rv I ry I 

iiv) — i{M) = / — — -— dx < „, , dx. 

^ ' ^ ' Jm f{x)+goi{x) -JmUx) 
Hence, i{y) converges to a finite limit as y — t- oo so that c explodes. 

□ 

2.5. Application of the analytic theory. We now pass to a probabilistic application of The- 
orem 1. 

Proof of Corollary 1. We consider first the case where Y is deterministic. Since Y is assumed 
to be strictly increasing, we can consider the unique non-negative stochastic process Z which 
satisfies 

Zt = x + Xj.t^^^^ + Yt. 

(The reader can be reassured by Lemma 5 regarding any qualms on measurability issues.) Since 
Z is non-negative. Theorems. 4.1 and 4.2 of Kallenberg (1992) imply the existence of a stochastic 
process X with the same law as X such that 

Zt = x+ [ Zl/^'dXs + Yt. 
Jo 

Hence Z is a weak solution to (6). 

Conversely, if Z is a solution to (6), we apply Theorems 4.1 and 4.2 of Kallenberg (1992) to 
deduce the existence of a stochastic process X with the same law as X such that 

Zt = x + Xjtz^^^ + Yt. 

Considering the mapping (/,<?) i— )• F(f,g) that associates to every admissible breadth-first pair 
the solution h to (7), we see that Z has the law of f(^X, x + Y^^ . Hence, weak uniqueness holds 
for (6). 

When Y is not deterministic but independent of X, we just reduce to the previous case by 
conditioning on Y (or by augmenting the filtration with the u-field criYt :t >{))). □ 
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3. CBI PROCESSES AS LAMPERTI TRANSFORMS 

We now move on to the analysis of Theorem 2. Let X and Y be independent Levy processes 
such that X is spectraUy positive and y is a subordinator under the probabihty measure P. Call 
^ and $ their Laplace exponents (taking care to have $ > as for subordinators). Note that the 
trajectories of Y are either zero, piecewise constant (in the compound Poisson case), or strictly 
increasing. 

Let Z be the stochastic process that solves 

Zt = x + Xjt^^^^ + Yt 

and has no spontaneous generation (when Y is compound Poisson). To prove that Z is a 
CBI(^, <1>) we should see that it is a cadlag and homogeneous Markov process and that there 
exist functions ut : (0, oo) — )■ (0, oo) and vt : (0, oo) — )■ (0, oo) satisfying 

.9^ i gutW = o utiX) (9_v,(X) = ^utiX)) 

^' \uo{X) = X \vo{X) = 

and such that for all A, i > 0: 

(At this point it should be clear that the equation for u characterizes it and that, actually, for 
fixed A > 0, 1 1— )• ut{X) is the inverse function to 

3.1. A characterization Lemma and a short proof of Lamperti's Theorem. The way to 

compute the Laplace transform of Z is by showing, with martingale arguments to be discussed 
promptly, that 

(10) E(|e-^^') = y E(^['f(A)Z, -$(A)]e-^^=) ds. 

We are then in a position to apply the following result. 

Lemma 4 (Characterization Lemma). If Z is a non-negative homogeneous Markov process with 
cadlag paths starting at x and satisfying (10) for all X > then Z is a CBI{^, <I>) that starts at 

X. 

Remark. Note that the hypotheses on the process Z of Lemma 4 do not allow us to use gener- 
ator arguments which would shorten the proof, for example by using the characterization of the 
infinitesimal generator of a CBI process through exponential functions. 

Proof of Lemma 4. Let us prove that the function 

G{s) = E(^e-"*-=(^)^--^*'=(^)) 

satisfies G'{s) = for s G [0,t], so that it is constant on [0,t] implying the equality 

E^e-^^*) = G{t) = G(0) = e-^"'(^)-^*W. 

We then see that Zt has the same one-dimensional distributions as a CBI(^', that starts at x, 
so that by the Markov property, Z is actually a CBI(^, $). 
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We now analyze both summands to later divide by h and let /i — )• 0. 
For the first summand, use (10) to get: 



G{s + h)- Ej^e"^''"'-^-''^^)-''*— 



-Vt-s-hW 



s+h 



E 



dr, 



so that, since Z has cadlag paths, we get 

lim i \g(s + /i) - Efe~^="*— '^(^)') 

jE(^g-n,_.(A)Z.-.,_.(A) [^^^ ^ _ ^ ^ ^,,_^(A)]) 



For the second summand in the right-hand side of (11), we differentiate under the expectation 
to obtain: 

lim —i^f g-'^t-s-hWZs-vt-s-hW _ g-Mt-s(A)Zs-ft-s(A) 
h^o h \ 



El g-"t-s{A)2s-<^t-s(A) 



■ dut-s{X) ^ dvt-s{X) 

Zjx 7^ r 



ds 



ds 



We conclude that G'{s) = for all s G [0,t] using (9). 



□ 



A simple case of our proof of Theorem 2 arises when Y = 0. Recall from Proposition 4 the 
notation 

Ct= Zgds. 
Jo 



Proof of Theorem 2 when $ = 0. This is exactly the setting of Lamperti's theorem stated by 
Lamperti (1967a). 

When $ = (or equivalently, Y is zero), then Ct is a stopping time for X (since the inverse of 
C can be obtained by integrating l/(x + X)). Since Z is the time-change of X using the inverse 
of an additive functional, Z is a homogeneous Markov process. (Another proof of the Markov 
property of Z, based on properties of IVP(X, x + Y) is given in (3) of Lemma 5.) Also, we can 
transform the martingale 

e-^-^t - ^(A) [ e-^^' ds 
Jo 

by optional sampling into the martingale 



-\Zt 



*(A) / e-^^'Zsds. 
Jo 



We then take expectations and apply Lemma 4. 



□ 
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3.2. The general case. For all other cases, we need the following measurability details. Consider 
the mapping Ft which takes a cadlag function / with nonnegative jumps and starting at zero, 
a cadlag g starting at zero (either piecewise constant or strictly increasing), and a nonnegative 
real x to c'^{t) where c solves IVP(/, x + g) and has no spontaneous generation (if g is piecewise 
constant). (Note that these conditions uniquely specify a solution to IVP(/, x + g).) Then 

(12) Zt+s = Ft{Xc,+. - Xc^ , - y„ Zs) . 

The mapping Ft is measurable. Indeed, we can view it as the composition of three measurable 
mappings. The first one is the mapping that associates to (/, g + x) the unique solution to 
IVP(/, 5) (without spontaneous generation), from the space of admissible breadth-first pairs 
equipped with the cr-fields generated by the projections (/, g) 1— )• f{t) and (/,<?) 1— )• g{t) for 
any t > to the space of non-decreasing continuous functions with cadlag derivative (equipped 
also with the a-field generated by projections). This mapping is measurable when g = Q hy 
measurability of the Lamperti transformation. Next, when g is piecewise constant this follows 
by concatenation of Lamperti transforms as in the introduction to Section 2, and for strictly 
increasing g this follows since the unique solution to IVP(/, g) is the limit of solutions to IVP(/, gn) 
with piecewise constant functions gn, as seen in the proof of Theorem 1. The second mapping 
sends a continuous function with cadlag derivative to its derivative, which is measurable by 
approximation of the derivative by a sequence of diff^erential quotients. The third mapping is 
simply the projection of a cadlag function to its value at time i; its measurability is proved in 
Theorem 12.5 p. 134 of Billingsley (1999). 

We suppose that our probability space (17, ^,P) is complete and let ^ stand for the sets in 
^ of probability zero. For fixed y, t e [0, 00], let = V V 

Lemma 5 (Measurability details). 

(1) The filtration (^y,y > O) satisfies the usual hypotheses. 

(2) Ct is a stopping time for the filtration (^y,y > O) and we can therefore define the a-field 

^^•, = {Ae^:An{Ct<y}e^;^}. 

(3) Z is a homogeneous Markov process with respect to the filtration (§f(^^,2/ > O) . 
Proof. 

(1) We just need to be careful to avoid one of the worst traps involving a-fields by using 
independence (cf. (Chaumont and Yor, 2003, Ex. 2.5, p. 29)). 

(2) We are reduced to verifying 

(13) {Ct <y}e ^l. 

We prove (13) in two steps, first when Y is piecewise constant, then when Y is strictly 
increasing. 

Let Y be piecewise constant, jumping at the stopping times Ti < r2 < . . ., and set 
To = 0. We first prove that 

(14) {C7T„<y}n{T„<t}el^,* 

and this result and a similar argument will yield (13). The membership in (14) is proved 
by induction using the fact that C can be written down as a Lamperti transform on 
each interval of constancy of Y . Let It be the functional on the subspace of Skorohod 
space consisting of functions with non- negative jumps that aids in defining the Lamperti 
transformation: when applied to a given function /, we first define 



ro(/)=inf{t>0:/(t)=0} 
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and then 

/■*ATo(/) 1 

Itif)= jT-^ds. 
Jo f{s) 

We then have 

{Ct, <y}n {Ti <t} = {ly{X + Yo) > Ti A t} n {Ti <t}€ 

If we suppose that 

{Cr„ < y} n {r„ < G 

then the decomposition 

oo 2-™L2'"(j/-9)J . , l._L n 

{CT„,,<y}n{r„+i<t}= U U U 2^<^^.<^ ni^"+i^*} 

ge(0,iy)nQm=l fe=0 ^ ^ 

allows us to obtain (14). Then, the decomposition 

oo oo 2-™L2'"{?/-9)J , i- , 1 ^ 

{c^<y^=\^ u u u {^"^^<wn^^^-.<^ 

gives (13) when Y is piecewise constant. 

When Y is strictly increasing, consider a sequence e„ decreasing strictly to zero and a 
decreasing sequence (vr^) of partitions of [0, t] whose norms tend to zero, with 

vr„ = {t^ = 0<t? <..-<t^„=t}. 

Consider the process y" = (XDsgiotj defined by 

^" = en + 5^ yt.i[,n_^_,„)(,) + yas=t. 

Since vr" is contained in vr"^^ and e„ > e„+i, > Y^^^. If C" is the solution to 
IVP(X, x + y") with no spontaneous generation (defined only on [0, t]), then Lemma 2 
gives C"' > C"'~^^. Hence, (C") converges as n — )• oo, and since the limit is easily seen to 
be a solution to IVP{X,x + Y); the limit must equal C by the uniqueness statement in 
Theorem 1. To obtain (13), we note that 

and 

{Ci<y} = U{Cr<y}. 

n 

(3) Mimicking the proof of the Strong Markov Property for Brownian motion (as in (Kallen- 
berg, 2002, Theorem 13.11)) and using (13), one proves that the process 

{Xct+s - Xct,Yt+s - Yt)s>o 
has the same law as {X, Y) and is independent of which we can restate as 
{Xct+s — Xct, YtJf-s — Yt)s>Q has the same law as {X, Y) and is independent of (X*^*, y*) 

where Xf* = Xc^as and Y^ = YtAs- 
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Equation (12) implies that the conditional law of Zt+g given 1^^^ is actually Zg measurable, 
implying the Markov property. The transition semigroup is homogeneous and in t units 
of time is given by the law Pt{x,-) of Ft{X,Y, x) under P. Note that this semigroup is 
conservative on [0,oo]. 

□ 

We will need Proposition 4 for our proof of Theorem 2. 
Proof of Proposition 4. Consider the filtration = (^t) given by 

= a{Xs : s < t) V a{Ys : s > 0) V =3^. 

If Y is strictly increasing, then C is strictly increasing and continuous. For fixed e > 0, let 
Ti < T2 < . . . be the jumps of Y of magnitude greater than e. Arguing as in Lemma 5, we see 
that C-p. is a ?^-stopping time which is the almost sure limit of the ?f -stopping times 
as n — >• 00. Since X is a 1^-Levy process and C|-^._^^^-)+ < Ct^ for all n, quasi left-continuity of 
X implies that X o C does not jump at Ti almost surely. Since this is true for any e > 0, then 
X o C and Y do not jump at the same time. 

If Y is compound Poisson, we argue on its constancy intervals, denoted [Tj_i, Tj), i = 1, 2, . . .. 
On the set {Cg < Ct^ for all s < Tj}, we can argue as above, using quasi left-continuity. On the 
set {Cg = for some s < Tj}, we note that X reaches —Yt^ for the first time at Ct^- The 
hitting time of — Ir^ by X is a 1^-stopping time which is the almost sure limit of the hitting times 
of — Iji + as n — )• 00. The latter are strictly smaller than the former since X has no negative 
jumps. Hence, by quasi left-continuity, X is almost surely continuous at Ct,. □ 



Proof of Theorem 2. Since 



y 

e-^^y - $(A) / e-^^' ds 



'0 / t>o 

is a (^^y) -martingale, it follows that M = (Mf)^>Q, given by 



Mt = e-^^c^t - ^'(A) / e'^^'^'^Zsds 
Jo 



is a (^^Q)^^Q-local martingale. With respect to the latter filtration, the stochastic process 
(^i)t>o given by 



Nt = e"^^* + $(A) f e-^^= ds 



(A) / 
Jo 



is a martingale. Hence e~^'^°'~" and e~^(^'+^) are semimartingales to which we may apply inte- 
gration by parts to get 

ft 



g-AXoC ^-Xx-XY 



e"^^* = Local Martingale + [ e'^^" [^(A) Zg - ^(A)] ds + 

Jo 

where the local martingale part is 

/* e-^(^'+^=) dM, + f e-^^°^^dNs. 
Jo Jo 

Since X o C and Y do not jump at the same time by Proposition 4 and Y is of finite variation, 
we see that 

0. 



g-AXoC ^-Xx-XY 



(Cf. (Kallenberg, 2002, Theorem 26.6. (vfi)). 
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We deduce that 



Jo 



is a martingale, since it is a local martingale whose sample paths are uniformly bounded on 
compacts thanks to the non-negativity of Z. Taking expectations, we get (10) and we conclude 
by applying Lemma 4 since Z is a Markov process thanks to Lemma 5. □ 

3.3. Translating a law of the iterated logarithm. 

Proof of Corollary 4. Let X be a spLp with Laplace exponent 4*, <1> be the right-continuous 
inverse of ^ and 

/,x log I log i I 

a(t) = . 

^{t-^ log |logt|) 

Recall that ^ is the Laplace exponent of the subordinator T = {Tt,t > 0) where 

Tt = mi{s>0:Xs< -t} . 

(Cf. (Bertoin, 1996, Ch. VII, Thm. 1).) If d is the drift coefficient of l> then Proposition 1 of 
(Bertoin, 1996, Ch. Ill) gives 

A— >oo A 

Hence, 

'a{t)^t/d ifJ>0 



as t 0+ j 
We now assert that if at —t- 1 as t —)■ then 



t = o{a{t)) ifd = 



t-s>o a{t) 

This is clear when d > 0, so suppose that d = 0. Since 1 1— t- log |logt| is slowly varying at zero, it 
suffices to show that if — ?■ 1 as A — >• oo then 

(15) hm = 1. 

^ A^oo $(A) 

However, concavity of increasingness and nonnegativity give (if 6 > 1) 

$(6A) /b < ^(A) < $(6A) 

which implies (15). 

As noted by Bertoin (1995), Fristedt and Pruitt (1971) prove the existence of a constant C 7^ 
such that 

liminf — = C. 
t^o a[t) 

Let Z be the unique solution to 

Zt = x + Xjt^^^^ + Yt 

with X > 0, where X and Y are independent Levy processes, with X spectrally positive of 
Laplace exponent ^ and Y a subordinator with Laplace exponent Since Zq = x, and Z is 
right-continuous, then 

1 /■* 

lim - / Z,,ds = X 
t^-o+ t 
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almost surely. Hence 



a [ Jq ds 

lim 



t^o+ a{xt) 
and so 



liminf y;^' =C 
t->o+ a(xt) 

On the other hand, if d is the drift of ^ then 

lim — = d, 

(cf. (Bertoin, 1996, Ch. Ill, Prop. 8)) so that if d = 0, Ft = o[a{t)) and 

liminf —. — — = C. 

i-!>o+ a{xt) 

If d > then by Proposition 8 by (Bertoin, 1996, Ch. Ill) we actually have 

liminf — = — 1 

t-s>o ait) 



so that 



liminf — % = — 1 + — . 

t^o+ a\xt) X 



□ 



3.4. Explosion criteria for CBI. As a probabilistic application of the deterministic explosion 
criteria of Proposition 3, we prove Corollary 5. 

Proof of Corollary 5. Let x > and consider a spectrally positive Levy process X with Laplace 
exponent ^ independent of a subordinator Y with Laplace exponent Let Z be the unique 
solution to 

Zt = x + Xjt^^^^ + Yt, 
which is a CB^^I', that starts at x. Also, let 

r-i 

Ct = I Z. ds. 





Let y be a non-zero subordinator. Path by path, we see that Z jumps to infinity if and 
only if either X jumps to infinity or Y does. However, the probability that either X or 
Y jumps to infinity is positive if and only if either ^'(O) > or <I>(0) > 0. When Y is 
zero, Z jumps to infinity if X jumps to infinity and never reaches — x, which has positive 
probability. 

The Ogura-Grey explosion criterion for continuous state branching processes (as stated 
just before Corollary 5) can be restated as follows: a CBI(^', 0) started at x reaches oo 
continuously at a finite time with positive probability if and only if ^(0) = and ^ is an 
explosive branching mechanism. It is also simple to see that a CBI(^',0) jumps to oo at 
a finite time with positive probability if and only if ^'(0) > 0. 

Path by path, we see that if Z reaches oo continuously (say at time r), then Y does not 
jump to infinity on [0, r) . Also, if we let C be the unique solution to IVP(x + Y^- + e + X,0) 
and Z as the right-hand derivative of C, where e > 0, then C < C on (0, r) (as follows 
from the argument proving Lemma 2). Hence C explodes on [0,r). We conclude that the 
branching mechanism of Z is explosive by the Ogura-Grey explosion criterion. Hence, 
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the assumption ¥{Z reaches oo continuously) > imphes that ^'(0) = and that ^ is an 
explosive branching mechanism. 

On the other hand, if ^(0) = and ^ is explosive, let $ = <I>(0), y be a subordinator 
independent of X with Laplace exponent so that sending Y to infinity at an exponential 
random variable with parameter ^"(0) (independent of both X and Y) leaves us with 
a subordinator with Laplace exponent <I> independent of X. Let be a solution to 
IVP(x/2 + X, 0) and be a solution to IVP(rE + X, Y) so that < and hence 
explodes if does. Let be the right-hand derivative of C*. is a CBI(\I',0) starting 
at x/2 while is a CBI(^', started at x; notice that the process Z obtained by sending 
Z'^ to infinity at the same exponential as Y leaves us with a CBI(^', <I>). By assumption, 
X cannot jump to infinity and Z^ explodes with positive probability. Hence, explodes 
with positive probability and can only do so continuously. Hence, 

F(Z reaches oo continuously) > e~*'^'^^P(Z^ reaches oo continuously by time i) 

and the right-hand side is positive for t large enough. 
(3) We also deduce that 



condition for the latter is that, additionally, ^ is not zero. Indeed, when $ is not zero 
then Ft — )• oo as t — )• oo. Since is explosive and $(0) = then limt_^oo Xf = oo and so 
Proposition 3 implies that the solution to IVP(X, x + Y) explodes. If $ = then Z'^ is 
a CBI{^, 0), which cannot explode continuously almost surely since the probability that 
Z^ is absorbed at zero is the probability that X goes below —x, which is positive. 



We now turn to the proof of Theorems 3 and 4, and of Corollaries 6 and 7, which summarize 
the stability theory for IVP(/, g). 

4.1. Proof of the analytic assertions. In order to compare the initial value problem IVP(/, g) 
with the functional inequality (5), we now construct an example of an admissible breadth-first 
pair {f,g) such that IVP(/, g) has an unique solution but (5) has at least two. Indeed, consider 
g = and take 



Then IVP(/, 5) has an unique solution, by Theorem 1, since / is strictly positive. The solution 
is the function c given by 



is also a solution to (5). Hence, the assumption of Theorem 3 is stronger than just uniqueness of 
IVP(/, 5) although related (as seen by comparing Theorem 1 and Proposition 1). 
We start with a proof of Proposition 1 . 




= 1. A necessary and sufficient 



□ 



4. Stability of the generalized Lamperti transformation 
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Proof of Proposition 1 . Let c be any non-decreasing solution to 

t f-t 
/_ o c(r) + g{r) dr < c(t) — c(s) < / / o c(r) + g{r) dr, for s < t 



such that c(0) = 0. This automatically implies continuity of c and so f o c + g is cadlag and does 
not jump downwards. 

Note that c is strictly increasing if /_ + g{0) is strictly positive or g is strictly increasing, we 
have equalities in (5), implying that c solves TVP{f,g) which has a unique solution with these 
hypotheses. Indeed, if /_ +^(0) is a positive function, then the lower bound integrand is strictly 
positive and so c cannot have a constancy interval. If on the other hand g is strictly increasing, 
note first that the non- decreasing character of c implies, through (5), that foc+g is non-negative 
(first almost everywhere, but then everywhere since it is cadlag). Also, foc + g can only reach 
zero continuously since it lacks negative jumps. If c had a constancy interval [s, t] with s < t, 
there would exist r S (s, t) such that 

/_ o c(s) + g{s) = f-o c(s) + g{r) = 

which implies that g has a constancy interval on [0, t], a contradiction. Hence, c has no constancy 
intervals. 

When g is a constant and f-+gis absorbed at zero then also f + g is absorbed at zero and at 
the same time. Hence, c is strictly increasing until it is absorbed, so that again both bounds for 
the increments of c are equal. Then c solves IVP(/, i^) which has an unique solution under this 
hypothesis. □ 

We now continue with a proof of Theorem 3. It is divided in two parts: convergence of 
the cumulative population which is then used to prove convergence of population profiles. The 
strategy is simple: we first use the functional equations satisfied by (cn) to prove that Cn /\ K is 
uniformly bounded and equicontinuous. Then, we pass to the limit in the functional equations 
satisfied by Cn to see that any subsequential limit oi Cn /\ K equals c A K. (This is where the 
assumption that 5 has an unique solution comes into play.) Having established convergence of 
Cn to c we then verify some technical hypotheses enabling us to apply results on continuity of 
composition and addition on adequate subspaces of Skorohod space and deduce that fn°Cn + 9n 
converges to f o c + g. 

Proof of Theorem 3, convergence of cumulative populations. Let K,e > and consider the se- 
quence Cn A K consisting of non-decreasing functions with cadlag right-hand derivatives. Since 



< D+Cn AK{t) = lc„{t)<K X 



[fn° Cni[t/an\(Tn) + gn{t)]^ if fT„ > 
fn°Cn{t)+gn{t) if fJn = 



< sup f{y) + g{t) + e 
y<K 

for n large enough (by the convergence oi fn ^ f on [0, K] with the Ji topology and gn ^ 9 
on [0, t] with the Ji topology), we see that the sequence c„ A K is uniformly bounded and 
equicontinuous on compacts. To prove convergence of c„ A -fT (uniformly on compact sets), it is 
enough to prove by Arzela-Ascoli that any subsequential limit is the same. Let t > and c„j. A K 
be a uniformly convergent subsequence on [0,t]. Denote by c its uniform limit, which is then 
non-decreasing. If s G [0,t] is such that c(s) < x then c„^(s) < x for A; large enough. Since / has 
no negative jumps then 

liminf/(x) = f-{y) and limsup/(x) = f{y) 
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SO that 

/_ o c < lim inf fn,^ o Cn^ < lim sup /^^^ o Cn,^ < / o c. 
By Fatou's lemma, for any si < S2 < s 

/ [f- o c{r) + g{r)]~^ dr < c{s2) - c{si) < / [f o c{r) + g{r)]^ dr. 

J Sl J Si 

As c is non-decreasing, we might remove the positive parts in the above display and conclude, 
from uniqueness to (5), that c = c on [0, s]. If, on the other hand, c(s) = K then c^j. AK{s) — )• K 
which implies that Cn^. A K ^ c A K uniformly on compact sets. 

Let r be the explosion time of c. If t < r then c(t) < oo and so (choosing K > c{t) in 
the paragraph above), we see that — )■ c uniformly on [0,t]. If t > r then c{t) = oo and so 
Cn{t) AK —f K ioi any K > Q. Hence c„(t) — ;> oo = c(t). □ 

Proof of Theorem 3, convergence of population profiles. Let = D^Cn and h = D^c. We now 
prove that /i„ — )• /i in the Skorohod Ji topology if the explosion time r is infinite. Recall that 
h = f o c + g and that 



hr. 



fn°Cn+ gn if (Jn = 

[fn o Cn{[t/an\crn) + 5n ( LV^nJ ^^n)] ^ if fT„ > 



Assume that cr„ = for all n, the case C7„ > being analogous. Then the assertion /i„ — t- /i is 
reduced to proving that: /n o Cn — >• / o c, which is related to the continuity of the composition 
mapping on (adequate subspaces of) Skorohod space, and then deducing that fn°Cn+9n f°c+g, 
which is related to continuity of addition on (adequate subspaces of) Skorohod space. Both 
continuity assertions require conditions to hold: the convergence fn ° Cn ^ f ° c can be deduced 
from (Wu, 2008, Thm. 1.2) if we prove that / is continuous at every point at which is 
discontinuous, and then the convergence of fn° Cn + 9n will hold because of (Whitt, 1980, Thm 
4.1) since we assumed that foe and g do not jump at the same time. Hence, the convergence 
hn ^ h is reduced to proving that / is continuous at discontinuities of c~^. 

If c is strictly increasing (which happens when g is strictly increasing or / + 5^(0) > 0), then 
is continuous. (This is the most important case in the stochastic setting, since otherwise 
immigration is compound Poisson, therefore piecewise constant, and one might argue by pasting 
together Lamperti transforms.) 

Suppose that c is not strictly increasing and let x be a discontinuity of c^^. Let 

s = c~^{x—) < c^^{x) = t, 

so that c = X on [s, t] while c < x on [0, s) and c > x on {t, oo). Since D+c = foc + g = Ooii 
[s,t), we see that g is constant on [s,t). We assert that 

inf {y > : f{y) = -g{s)} = x. 

Indeed, if / reached —g{s) at x' < x, there would exist s' < s such that 

= / o c(s') + gis) > f o c{s') +g{s') > 

so that actually g is constant on [s',t). Hence, c has spontaneous generation which implies there 
are at least two solutions to IVP(/, 5): one that is constant on {s',t], and c. This contradicts the 
assumed uniqueness to (5). Since / has no negative jumps and reaches the level —g{s) at time x 
then / is continuous at x. 

Finally, we assume that the explosion time r is finite but that /n(x) ,/(x) — 00 as x — t- 00 
uniformly in n and prove that — ?■ /i in the uniform Ji topology. Let e > 0, d be a bounded 
metric on [0, 00] that makes it homeomorphic to [0, 1], and consider M > such that d{x, y) < e 
if X, y > M. Let K > he such that /(x) , fnix) > M if x > K and n is large enough. Let T < t 
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be such that / is continuous at c(T) and K < c(T). Then, fn^fin the usual Ji topology on 
[0, c(r)] and, arguing as in the non-explosive case, we see that 

K = fn ° Cn + Qn ^ f ° C + g = h 

in the usual Ji topology on [0,T]. Hence, there exists a sequence (An) of increasing homeomor- 
phisms of [0, r] into itself such that — h o — )• uniformly on [0,T]. Define now A„ to 
equal A„ on [0,r] and the identity on [T, oo). Then (A^) is a sequence of homeomorphisms of 
[0, oo) into itself which converges uniformly to the identity, and since K < c{T), then K < Cn(T) 
eventually and so M < h eventually thanks to the choice of K, so that d{hn{t) , h{t)) < e on 
[T, oo) eventually. Hence, /i„ — )• /i in the uniform Ji topology. □ 

In order to apply Theorem 3 to Galton- Watson type processes, we need a lemma relating the 
discretization of the Lamperti transformation and scaling. Define the scaling operators by 

S'afit) = If (at) . 

Let also c*^ be the approximation of span a to YVF{f,g), which is the unique function satisfying 

c''{t)= [\foc^als/a\)+g{als/a\)]+ ds. 
Jo 

We shall denote c'^{f, g) to make the dependence on / and g explicit in the following lemma and 
denote by h'^{f,g) the right hand derivative of c'^{f,g). 

Lemma 6. We have: 

Syif, g) = c-l'^ [s'/y, S'J'^g) and S'J'^h'^if, g) = h"''^ (^s'/^, S'J'^g) . 

The proof is an elementary change of variables. 
4.2. Weak continuity of CBI laws. 

Proof of Corollary 6. Let Xn and X be spLps with Laplace exponents and ^ and and Y 
be subordinators with Laplace exponents and $ such that Xn (resp. X) is independent of Yn 
(resp. Y). 

The hypotheses — ^ ^ and — )• <I> imply that {Xn,Yn) converges weakly to {X,Y) in the 
Skorohod Ji topology. By Skorohod's representation theorem, we can assume that the conver- 
gence takes place almost surely on an adequate probability space. 

Let Zn (resp. Z) be the Lamperti transform of {Xn,Xn + Yn) (resp. {X,x + Y)). When X is 
non-explosive, Propositions 2 and 4 and Theorem 3 then imply that Z„ converges almost surely 
to Z, which is a CBI(^', $) thanks to Theorem 2. 

When X is explosive, let p be a distance on [0, oo] which makes it homeomorphic to [0, 1] and, 
for any e > 0, choose Me such that p{x,y) < e if x,y > M^. Recall that doo stands for the 
uniform Ji topology. Since the X^ — )• X and y" — )• y in the usual Skorohod topology as n — )• oo 
almost surely, then reasoning as in the proof of uniform Ji convergence of Theorem 3, we see 
that, for any e > 0, 

P(doo(^", Z) > e, X^, Xs > Me for ah s > f) ^ as n ^ oo. 

However, choosing t and M big enough, we can make 

P(X," < M for some s>t) 

arbitrarily small for all n large enough, so that doo{Z", Z) — )■ in probability, which is enough to 
guarantee that — t- Z weakly in the uniform Ji topology. Indeed, since X is explosive, we have 
that ^-'(0+) = -oo (cf. (Lambert, 2008, Proof of Theorem 2.2.3.2, p. 95)) which means that X 
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drifts to oo (cf. Corollary 2.ii in (Bertoin, 1996, Chapter VII)). Since the latter result implies 
that the negative of the infimum of X has an exponential distribution of parameter ry, where 

r? = inf {A > : ^(A) = 0} , 

we see that 

F{Xs < M for some s > t) < F{Xt < 2M) + F{Xt > 2M and Xs < M for some s > t) 

< F{Xt < 2M) + 6"''*^ 

Since X drifts to infinity, the term F{Xt < 2M) goes to zero as t — )• oo. Asymptotically, the same 
bounds hold for X" since — )• ^ and hence, by convexity of ^, 

lim (inf {A > : ^'"(A) = 0}) = inf {A > : ^(A) = 0} = r/. 
n— >oo 

□ 

4.3. A limit theorem for Galton- Watson processes with immigration. 

Proof of Corollary 7. By Skorohod's theorem, if X and Y are Levy processes whose distributions 
at time 1 are ji and v then: 

^ X and 5." Y, 

where the convergence is in the Ji topology. Assume first that X is non-explosive. 
We can apply Lemma 6 to get either 



or 



Let Z be the unique process satisfying 
as in Proposition 2. If ^^\knlx\l (^\knlx\ — c G [0,oo), we see that 
thanks to Propositions 2 and 4 and Theorems 2 and 3. 



When ^^[fc„/a;j/o[fc„/a;J ~^ oo, let Z instead be the unique solution to 



Zt = x + Xjt^^^^. 



Then 

c 



5^;"/^ ^ Z. 



When X is explosive, the arguments in the proof of Corollary 6 show that, in order to obtain 
the stated convergence in the uniform Ji topology, it is enough to prove that for all M > 

lim lim limsupPf — i < M for some s > t | =0. 

Since X drifts to infinity if it is explosive, ^ has an unique positive root which we denote rj. 
Let 

G„(A)=Efe-^^r). 
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Recall that since the increments of X"^ are bounded below by —1, minus the random variable 



In = minX^ 

m>0 



has a geometric distribution with parameter e"''" where rjn is the greatest non-negative real num- 
ber at which Gn achieves the value 1 (cf. Corollary 5.5 in (Asmussen, 2003, Part B, Ch. VIII§5, 
p. 235) or the forthcoming Lemma 7). By log-convexity of Gn-, rjn = inf {A > : Gn[X) > 1}. If 
we assume the convergence of nr]n to 77 as n — )• c«, we see that 

limsuppf > M] = e-''^^ 



n 



We now use the Markov property to conclude that if the distribution of Xi is continuous at M, 
then 

limsuppf < M for some s > t] < ¥{Xt < 2M) + ¥{Xt > 2M) e'"'^^ . 



To conclude, we should prove that nr]n — )• ij. This however is implied by the following conver- 
gence of Laplace transforms: 



Indeed, recall that E(e ^^^^ < 1 exactly on (0, t]) and that E^e '^/''^^"nj < 1 exactly on (0, nr/„). 

If we consider X < rj then E^e"^/""''"""^ < 1 for large enough n, so that A < nry„ for large enough 

n. This implies r] < liminfnn?7„; the upper bound is proved similarly. Convergence of Laplace 
transforms is actually the condition imposed by Li (2006) to prove limit theorems for Galton- 
Watson processes with immigration. That this already follows from our hypotheses is the content 
of the following lemma, which concludes the proof of Corollary 7. □ 

Lemma 7. Let X^ be a sequence of random walks with jumps in {—1,0,1,...} satisfying the 
conditions of Corollary 7 and suppose that X is not a subordinator. Then 

for all A > 0. 

This is the content of Theorem 2.1 of Grimvall (1974), we present a proof using basic fluctuation 
theory for independent increment processes. 

Proof. Using Skorohod's theorem again, we assume that j/n converges almost surely to X 
in the Skorohod Ji topology. Also, enlarge the probability space so that it admits an exponential 
random variable Rx of parameter A which is independent of X and X" . 
Let 

G„(A)=E(e-^^r), 

Since X is not a subordinator, then ¥{Xf = — 1) > for large enough n and we can assume that 
this happens for every n. Hence, G„(A) — )• 00 as A — )• 00 and we can define 

Fn{s) = inf {A > : G„(A) > l/s} for s G (0, 1]. 

Using optional sampling at the first time TJ^ at which reaches —k for the first time, applied 
to the martingale 

e-'^'^GniX)-"' , 

we obtain 



E{s^l^) 
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for s G (0, 1]. Define the random variables 

Ix = min Xs and I\ = min —XV^ ^< 
s<Rx s<Rx n LQ-n'^J 

Since [anR\\ has a geometric distribution of parameter e~^l°"^ , it fohows that 

so that —nl^ has a geometric distribution. Also, from Corollary 2 in (Bertoin, 1996, Chapter 
VII), Ix has an exponential distribution of parameter $(A) where 



$(A) = inf |a > : ^(^a) > a} . 



However, since X does not jump almost surely at Rx and the minimum is a continuous functional 
on Skorohod space (on the interval [0,/?^])) we see that /" converges weakly to Ix- This implies 

nF„(e-^/"") ^ $(A) 

and by passing to inverses, we get 

G„(A/n)°" ^ e*(^) 

for A > l>(0). 

Finally, if A G (0, ^(0)], pick p > 1 such that pX > ^(0); we have just proved that the sequence 

GnipX/n)"" n > 1, 
being convergent, is bounded. Hence the sequence 

e-^/"^-, n>l 
is bounded in Lp and converges weakly to e~'^^^. We then get 

□ 

4.4. A limit theorem for conditioned Galton- Watson processes. 

Proof of Theorem 4. Let be a Galton- Watson process with critical offspring law // such that 
= kn and is conditioned on Yli^i ^? — ^- Then, has the law of the discrete Lamperti 
transformation of the n steps of a random walk with jump distribution fi (the shifted reproduction 
law) which starts at and is conditioned to reach —kn in n steps; call the latter process X^, so 
that 

Z"" = h\kn + X'',0). 



Thanks to Chaumont and Pardo (2009), if kn/an — I then 
Thanks to Lemma 6, we see that 

Z" = /i""/"(5""X",0). 

Let a G (1, 2] be the index of the stable process in the statement of Theorem 4 and recall that, 
an is of the form n^l°'L{n) where L is a slowly varying function, so that a„ = o(n). Since is 
absorbed at zero (as easily seen by the pathwise construction of F^ by (Chaumont and Pardo, 
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2009, Thm. 4.3)), then Proposition 1 guarantees that the Lamperti transform Z of is the 
unique process which satisfies: 

^/J Zudu- — Zr dr < Fp^ Zudu' 



Theorem 3 imphes that 

n/an 

□ 
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